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A LOW  COST  NAVIGATION  MICROPROCESSOR  SYSTEM 


S.  Suraratrungsl  and  P.  S.  Noe 


ABSTRACT 

The  Global  Positioning  System  (GPS)  is  a satel- 
lite-based radio  navigation  system  currently  under  de- 
velopment by  the  Department  of  Defense.  This  paper 
describes  the  feasibility  of  a microprocessor  implemen- 
tation as  part  of  a low-cost  GPS  receiver.  A noisy 
model  is  considered  in  the  dynamic  system  simulation. 

A total  of  726  fixes  are  taken  during  a 6-hour  flight 
with  an  RMS  error  of  31.59  ft.  Results  are  compared 
with  a similar  run  on  an  Amdahl  470V/6. 

INTRODUCTION 


of  i satellite  (as  adetermined  from 
the  ephemeris  data) 

u.  = three  receiver  coordinates  (as  deter- 
J mined  from  the  estimated  latitude,  lon- 
gitude, and  altitude  of  the  receiver) 

r^  ■ range  from  the  iC^  satellite  to  the 
user 

b ■ clock  bias 

Rearranging  (1)  and  solving  for  r^  gives 


ri ' [(xirui 


)2  + (xi2~u2^2  + ^ 


.2 

xi3'u3) J 


The  NAVSTAR/Global  Positioning  System  (GPS)  is  a 
satellite-based  radio  navigation  system  authorized  for 
development  by  the  Department  of  Defense  in  December 
1973.  Early  in  the  1980s  a complete  system  consisting 
of  24  satellites  will  orbit  the  earth  in  circular  12- 
hour  orbits  at  an  altitude  of  11,000  nautical-miles  and 
inclined  at  63°  to  the  equator.  The  satellites  will 
transmit  pseudo-random  noise  codes  and  ephemerides  to 
the  users  worldwide.  The  user,  equipped  with  a small 
receiver,  receives  and  processes  the  navigation  signals 
from  which  he  will  be  abie  to  determine  his  location, 
velocity,  and  system  time. 

Hard-wired  processors  and  minicomputers  have  been 
used  effectively  as  the  major  processing  tool  for  sat- 
ellite navigat  Lon  in  the  past,  but  high  cost  and  weight 
impose  a practical  limitation  for  these  systems  to  cer- 
tain classes  of  users.  Recently  the  advent  of  low- 
cost,  high-performance  large-scale  integrated  micro- 
processors have  eliminated  this  limitation.  Microproc- 
essors have  several  advantages  over  hard-wired  proces- 
sors such  as:  low  cost,  light  weight,  less  power  con- 

sumption, reliability,  and  flexibility.  These  advan- 
tages make  it  feasible  to  develop  a microprocessor- 
based  design  for  a low-cost  GPS  receiver.  It  is  the 
objective  of  this  paper  to  show  the  feasibility  of  such 
a GPS  receiver. 

THE  ALGORITHM 

The  fundamental  problem  in  the  implementation  of  a 
GPS  low-cost  receiver  is  described  in  [21.  This  prob- 
lem involves  the  computation  of  user  position  compo- 
nents u. , U£»  u^,  and  clock  bias  b which  can  be  ob- 
tained rrom  the ^spherical  navigation  equations 


Obtaining  a Taylor  series  expansion  gives 
l I 2 

r.  = i.  + (3r./3U)  + (32r./3U2)  £ U+  ....  (3) 

i i i U i |U 


where  U is  the  user's  estimate  of  U and  r.  is  computed 
by  substituting  U into  (2). 

Linearizing  (3)  by  retaining  only  the  linear  terms  gives 
Sr.  - [dr./du,  dTj/du^  dr.fdu  3r./3b]j~U  (4) 

1 ' 1 A A A 1 J 1 , L' 

where  6r . = r . - r and  <5U  * U - U 

1 1 (5) 

Now  define  r^  * h^6U 

where  h^  is  the  row  vector  in  (4)  given  by 

ht  * jla1-x11)/(ri-b)  (u2-xi2)/(r1-b)/(u3-.xl3)/ 

(Vb)  1]  (6) 


5R  - 2 <5u  - HSU 

h. 


where  H is  the  4x4  matrix  of  partlals  of  R with  respi  ■ 
to  U,  evaluated  on  U,  and  R =[r,  r,  r»  r |T. 


The  desired  position  correction  are 


where  x^ , J-l,  2,  3 • three  known  position  components 


SU  - H SR 


The  update-estimated  position  is 
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Convergence  is  controlled  and  determined  bv  conparln 
6U  with  an  arbitrary  minimum  value,  e. 

The  Hotelling  algorithm  [3]  was  proposed  In  [2]  to  ,M 

used  to  Iteratively  update  H 1 In  (8)  as  the  nujof  1 
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In  solving  for  successive  GPS  position  fixes  [2].  This 
feature  Is  Implemented  as  step  (6)  below. 

In  summary,  the  algorithm  suggested  In  [2]  proceeds  as 
follows: 

1)  Guess  Initial  U;  select  A satellites;  compute 


2) 

3) 

A) 

5) 

6) 


G-H_1. 
Obtain  r 
i«l-*A.  _ 

Calculate  r 


1 and  xij 


from  receiver  data  for 
for  all  i. 


Obtain  6R  ■ R-R. 

Obtain  H matrix  from  (6). 


f 2I-HG  , 1 

If  G-G  . 

L J 

m m-1 

repeat  step  6);  otherwise,  go  to  7).  Here 

G . is  the  last  estimate  of  the  inverse  of  H 
m— l 

and  I is  the  identity  matrix. 

7)  Obtain  6U*G  5R. 

m 

8)  Update  position  estimate  U^U^  Repeat 

step  3)  through  8)  if  6U>t2»  a specified  tol- 
erance; otherwise  return  to  step  2). 


SIMULATION  MODEL 

Satellite  position  and  mission  profile  data  for  a 
USAF  C5A  aircraft  are  generated  in  realistic  computer 
simulation  programs,  SATGEN  and  PROFGEN,  respectively. 
The  programs  were  developed  at  the  Air  Force  Avionics 
Laboratory,  WPAFB,  Ohio. 

Satellite  Position  Generator  (SATGEN) 

The  Satellite  Position  Generator  (SATGEN)  simu- 
lates the  propagation  of  24  GPS  satellites  in  3 equally 
spaced  orbit  planes  at  63°  inclination.  Each  orbit 
plane  contains  a ring  of  8 satellites  which  are  equally 
spaced  45°  apart.  The  satellites  in  each  of  3 orbit 
planes  are  also  positioned  so  that  the  satellites 
crossing  the  equatorial  plane  (ascending  nodal  cross- 
ing) are  15°  apart. 

The  SATGEN  program  contains  an  algorithm  that 
tests  which  satellites  are  in  view  with  respect  to  the 
user’s  estimated  position,  and  subsequently  four  satel- 
lites are  selected  for  navigation  based  on  the  estab- 
lished criterion,  minimum  Geometric  Dilution  of  Pre- 
cision (GDOP)  [4J.  Currently  two  satellite  selection 
routines  used  in  [2]  were  considered-OPTSEL  and  SATSEL. 
The  OPTSEL  routine  computes  the  GDOP  for  all  possible 
combinations  of  satellites  that  are  in  view  and  selects 
4 satellites  with  minimum  GDOP.  The  SATSEL,  on  the 
other  hand,  selects  three  satellites  among  those  satel- 
lites in  view.  The  satellites  selected  are  those  each 
of  which  has  a largest  component  along  x,  y,  and  z co- 
ordinates respectively  (earth  centered,  earth  fixed  co- 
ordinate system).  The  fourth  satellite  is  chosen  such 
that  the  combination  with  previously  chosen  three  satel- 
lites will  result  in  minimum  GDOP.  Although  the  SATSEL 
routine  is  a suhoptimal  procedure,  it  provides  computa- 
tional advantages  over  the  optimal  GDOP  routine,  OPTSEL. 
For  example,  if  n satellites  are  in  view, |n^  GDOP  compu- 
tations are  required  to  determine  the  optimal  set  of 
satellites.  On  the  other  hand  the  suboptimal  method  re- 
quires only  n-3  GDOP  computations.  The  SATSEL  routine 
was  used  in  the  simulation. 


Mission  Profile  Generator  (PROFGEN) 

The  C5A  mission  profile  data  representing  a simu- 
lated 6-hour  flight  from  Travis  AFB,  California,  to 
Hickman  Field,  Hawaii,  were  generated  by  the  USAF 
Avionics  Laboratory,  WPAFB,  Ohio.  Those  data  arc  ref- 
erenced by  two  angular  coordinates  (latitude-longitude) 


and  the  altitude  above  the  reference  ellipsoid  (Inti— 
tude-longitude-alt itude  coordinate  system).  Since  cal- 
culations in  space  are  best  carried  out  in  earth  cen- 
tered, earth  fixed  (ECEF)  coordinate  system,  the  USAF 
C5A  mission  profile  data  were  transformed  to  ECEF  co- 
ordinate system  by  the  PROFGEN  routine. 

SIMULATION  RESULTS 

The  position  fixing  algorithm  [2]  is  implemented 
on  the  Intel lec  8/MOD  80  microcomputer  development  sys- 
tem which  utilizes  an  Intel  8080  microprocessor  with  an 
instruction  cycle  time  of  2 usee.  The  algorithm  re- 
quires 7k  bytes  of  read-only  memory  (ROM)  for  program 
storage  and  4k  bytes  of  random  access  memory  (RAM)  for 
data  storage.  A USAF  C5A  flight  is  simulated  with 
noise  in  the  data  generated  by  the  satellite  and  mis- 
sion profile  routines.  Position  fixes  are  obtained 
every  30  seconds  and  the  satellites  are  selected  every 
15  minutes.  A total  of  726  fixes  are  taken  during  a 6- 
hour  flight  with  ft.  in  step  8)  of  the  algorithm 

as  the  convergence  criterion.  With  this  value  of  e7 
only  one  iteration  is  required  for  each  fix,  and  the 
errors  are  220  ft.  or  less  as  shown  in  Figure  1. 


A similar  run  on  the  Amdahl  470V/6  computer  produces 
comparable  results  as  shown  in  Figure  2.  Duplicate  re- 
sults are  not  achieved  due  to  different  pseudo-random 
number  generator  routines  used  in  the  noise  model  for 
the  two  machines. 


It  should  be  pointed  out  that  the  Amdahl  470V/6  com- 
puter, used  in  the  simulation  along  with  the  Intel lec 
8,  is  for  comparison  purposes  onlv. 


■ 
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Computation  time  on  the  Intellec  Is  10  s/fix.  Aa 
far  as  the  speed  is  concerned,  fixing  time  on  the  In- 
tellec Is  quite  slow  but  perhaps  suitable  for  low-cost 
GPS  applications.  The  8-blt  microprocessor  is  a handi- 
cap for  this  system  if  the  10  s/fix  rate  is  considered 
to  be  too  slow. 


CONCLUSION 

It  has  been  shown  that  the  simulation  results  from 
the  microprocessor  system  are  comparable  with  those 
from  a general  purpose  computer  as  far  as  the  accuracy 
is  concerned.  Processing  speed,  on  the  other  hand,  is 
quite  slow  on  the  8-blt  machine  but  can  be  improved. 

The  following  are  some  techniques  which  are  recommended 
for  improving  the  processing  power  of  the  system: 

I.  A 16-bit  microprocessor  can  be  used  instead  of 
an  8-bit  machine  for  this  system.  Because  its  16-bit 
capability  can  process  16  bits  at  a time,  it  can  pro- 
vide greater  speed  than  the  system  built  with  an  8-bit 
microprocessor.  Moreover,  a 16-bit  system  can  work 
with  shorter  programs  using  less  memory. 

•2.  Since  most  microprocessor  speed  is  usually 
faster  than  memory  access  time,  the  processing  speed  is 
mostly  limited  by  the  speed  of  the  memory.  To  compen- 
sate for  this  mismatch  a high-speed  memory  can  be  con- 
sidered. 

3.  Use  the  hardware  multiply/divide  unit  to  proc- 
ess the  data  Instead  of  using  the  shift  and  add  algo- 
rithm from  the  software.  This  saves  programming  steps 
and  hence  Increases  processing  speed. 

A.  Utilize  multiprocessor  structure  in  the  sys- 
tem where  each  independent  task  is  handled  by  a dedi- 
cated microprocessor  or  set  of  microprocessors. 

5.  Alternatively,  the  bipolar  bit-slice  micro- 
processor can  offer  processing  power  which  is  far 
greater  than  that  presently  available  from  MOS  micro- 
processors. By  packing  their  processing  power  on  sev- 
eral matched  LSI  chips,  they  are  easily  expandable  to 
16-bit  or  even  32-bit  word  lengths,  and  they  can  be 
microprogrammed  to  handle  the  most  powerful  high-level 
Instruction  set  available.  However,  a bit-slice  micro- 
processor is  considered  more  expensive  than  those  of 
MOS  microprocessors.  A trade-off  between  speed  and 
cost  should  be  carefully  considered  in  the  design. 

The  discussion  and  recommendations  outlined  above 
leads  to  the  conclusion  that  a microprocessor  or  a set 
of  microprocessors  is  a feasible  way  to  implement  the 
processing  part  of  a low-cost  GPS  receiver. 
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